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1. INTRODUCTION 

Although the literature concerning motorcycle and bicycle 
dynamics is not large, such research usually has been concerned 
with the uncontrolled single-track vehicle* The work of Sharp 
[ij represents the most complete theoretical analysis of the 
uncontrolled motorcycle presently available, and serves as a 
theoretical basis for the study described herein. Sharp's analysis 
included roll, yaw, lateral translation, and steering degrees of 
freedom, and the resulting equations are linear with constant 
coefficients. For use here, tire aligning moments due to tire 
sideslip have been added to Sharp's equations. 

Of the few existing studies of a man and single-track vehicle 
as a closed-loop control system, the most significant are 
Van Lunteren, et al. (e.g., [2, 3]) and Weir [4]. While the 
former research represents a pioneering effort in obtaining rider 
transfer functions experimentally (with a bicycle simulator), it 
should be pointed out that Van Lunteren's major interest was the 
performance of the human operator under various conditions (drugs, 
etc.) and not the dynamics of the bicycle. Thus, the accuracy of 
the simulator dynamics with respect to real bicycles and the 
validity of the assumption of steering angle control (rather than 
steering torque) are questionable. On the other hand, the work of 
Weir, although strictly theoretical, was based upon the motion 
equations and motorcycle parametric data of Sharp, and represents 
a comprehensive study of the man-motorcycle system. 

It is the purpose of the work described here to present 
experimental results and relate them to the theoretical studies 
of Sharp and Weir. 


2, DESCRIPTION OF EXPERIMENTS 

The particular rider task under consideration here is stabili- 
zation of the vehicle roll angle by means of steering torques 
applied to the handlebars. Path following and control by body 
movements are not being studied. From theoretical considerations, 
Weir [4] has shown that roll stabilization is best accomplished by 
steering torque control and that body lean control is more suited 
to the path following task. An "inner loop" of the multi-loop 
man/motorcycle control system postulated by Weir is shown in 
Figure 1, in which 4 is the roll angle, t^ is the rider-applied 
steering torque, i includes wind and road disturbances, and n is 
the rider's remnant. The block diagram of Figure 1 represents the 
man-motorcycle system model that is used in the present study. 
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Figure 1. Block diagram of the man-motorcycle system* 



The motorcycle transfer function, or controlled element ,Yc(ju), 
can be calculated from the equations of motion of Sharp (!]. It 
will be shown that boi.h and the rider transfer function, 

Yp(jw), can also be obtained experimentally. 

To restrict the experiments to the roll stabilization task, 
road tests were performed on a 0.5-mile section of essentially 
straight paved road in good condition. The test motorcycle was 
operated next to an automobile, which carried the power supplies 
and recording equipment. Figure 2 shows the experimental situation, 
including the rigid brace mounted on the motorcycle to reduce rider 
body lean, and the third wheel that was used to measure roll angle. 
(Also visible is a horizontal arm supporting two weights which was 
removed for these experiments.) Additional instrumentation 
included a rate gyro to sense roll rate, a steering angle potentio- 
meter, and a strain-gaged bar to record steering torques. With 
this torque bar, it was necessary for the rider to steer with one 
hand. Throttle control was relocated to the rear frame. The 
effect of this unusual control arrangement upon the generality of 
the results is not known, hut it is felt that any such effects were 
small and, at least, considerably smaller than they would have 
been if the rider had been using position control instead of torque 
control . 

Three experienced riders, denoted Rider A, Rider B, and Rider 
C, were tested. These riders were chosen to be of approximately 
the same weight (170-196 lb.) to minimize changes in Yc(jw) due 
to variations in rider mass and mass distribution. 

A total of fifteen experiments were performed with Rider A 
{the author). These experiments took place on two days about four 
months apart. On the first day, experiments began with a 30 mph 
test, followed by a 15 mph test, followed by a 30 mph test, etc,, 
until eleven total were recorded, with three of each speed being 
suitable for analysis. The remaining nine tests wer^ selected 
from a group of ten (all at 30 mph) performed on the second day. 

On a single day, a total of thirteen tests were carried out 
with Rider B, with twelve of these being suitable for analysis. 
Similarly, on a single day, twelve tests were recorded with Rider 
C, eleven tests being suitable for analysis. Both riders were 
instructed to ride at a constant 30 mph while minimizing body 
movements with the aid of the brace attached to the motorcycle, 

A few practice runs were made before recording was begun. 

The duration of eacn 30 mph test was about 50-60 seconds, 
while the 1» mph tests lasted about 65 seconds. 

A major objective of the experiments was to identify the rider 
transfer function, Yp(jw). Since the wind/road disturbance i is 
unknown, it is necessary to estimate Y (jw) by an open-loop method. 


To reduce bias errors arising from a correlation between n and 
a method developed by Wingrove and Edwards [5, 6 , 7, 8 ] must be 
employed. When the remnant autocorrelation function is sufficently smal' 
for lags less than or equal to the rider’s time delay, ip, this 
method can remove most of the bias errors by first identifying 

e^^**^Yp(jw) , the transfer function having -♦(t-X) as input and 
t(S(t)*^as output, where 0 < X i Tp, since there is much less 
correlation between 0 (t-X) and nXt) than between ^ft) and n(t). 

The Wingrove-Edwards method must be coupled with an identification 
procedure that is constrained to identify only physically realizable 
systems. The impulse response method [ 6 , 9], which first identifies 
a discrete form of the impulse response function of the unknown 
system, and then determines the transfer function by Fourier 
transformation, is such a procedure and has been used in this study. 

3. DISCUSSION OF RESULTS 

Preliminary analysis of the experimental data for 30 mph tests 
indicated that the theoretical controlled element, Y(.(Jw), could 
be accurately identified with the following relationship: 






Where Y^(jo)) is the estimate of Y^(Ja>), is the estimated 

i ross- spectrum between ^ and t^, and estimated 

power spectrum of t^. An example of such an identification is 
shown in Figure 3. The theoretical Y^(jw), calculated from the 
eouations of motion and the motorcycle parametric data measured 
foi the test vehicle (Appendix A), is shown for comparison. The 
identification is seen to be the closest to the theoretical con- 
trolled element for frequencies of about 1-4 radians/second, where 
most of the power of the roll angle and steering torque records is 
located. 

The fact that the theoretical Yr(jw) can be identified experi- 
mentally provides support for the motorcycle equations of motion 
and also indicates that the rider's remnant, n, is a much greater 
source of excitation to the man-machine system than the wind/road 
disturbance, i. When the remnant is the primary disturbance to 
the system, the Wingrove-Edwards method is required to identify 
Yp(jw). 
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Figure 3. Identification of 30 »pJ>* 


Since the inpulse response and Wingrove-Edwards methods are 
not nearly as well documented in the literature as spectral analysis 
techniques, some rules regarding the interpretation of Bode diagrams 
of Yp(ju) were derived from applying the methods to the identifi- 
catiSn of known systems from artificially created data,^as outlined 
in Appendix B, A further aid to the inteiiiretation of Y^(jw) 

Bode diagrams was obtained by assuming; the remnant power^spectrum 
to be "white** and calculating the expected power spectrum of tx, 
based on the estimate of YpCJw) and the theoretical Yr(jw) 

(including the effect of the analog filter applied to the data] • 

Two example identifications of Yp(Ju) are shown in Figure 4. 

In general, the form of Yp(jw} for the 30 mph tests was found to 
be ^ 

■ • KpO (TjjJo) ♦ 1). (2) 

The value of Tji was negligibly small, except for the tests 
carried out with Rider A on the second day* Thus, for most of 
the 30 mph tests, a rider transfer function of the form, 

Yp(jw) = - Kp e P , (3) 

fits the data* In the three IS mph tests the rider transfer 
function was of the form 

" * Kp e ^P^“ • 1). (4) 


Estimated mean values of the parameters in Equations (2-4) are 
presented in Table 1* The confidence intervals for the true mean 
values were calculated on the assumption that the estimates of 
the means were unbiased* Since this assumption may not be valid, 
more certainty may be attached to the width of the confidence 
interval than its midpoint* 

The resulting rider transfer functions were found to be 
consistent with the crossover model, which was employed by Weir 
[4] in his theoretical analysis of the man*motorcycle system* The 
crossover model requires that 
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Fl^rare 4a. Identification of Y (Jw) for a 
single test (Rider 30 raph, 
first day). 
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Figure 4b. Identification of Yp(J«) for a 

single test (rider A, 30 n^h, 
second day; the case In which 
i evidenced the most lead). 
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TABLE 1 


SUMMARY OP RIDER TRANSFER PONCTIONS 



Day of Test 

1 

1 

2 

3 

i| 


Speed, mph 

15 

30 

30 

30 

30 


Rider 

A 

A 

A 

B 

C 





— — 4 • • 


-T-jw 


Form of 

-Kpjue P 

-Kpe P 

-V 


-Kpe P 


YpUo.) 

•(TjjJw + 1) 


•(TjjJw + 1) 


hean Values 
pf Estlmi.ted 

^p* radian 

75.6 

277 

26l 

167 

206 

{Quantities 

Tp, seconds 

0.300 

0.398 

0.300 

0.328 

0.299 


, radians 
secon<r 

9.0J4 


6.«I5 




degrees 

25 

38.3 

ii8.6 

41.4 

39.1 

90% 

Confidence 

« Ib-ln 

*p* radian 

(53.6, 97.6) 

(219, 335) 

(236, 286) 

(152, 183) 

(185, 227) 

Intervals 

seconds 

P 

(.316, .38it) 

(.360, ,si6) 


(.387, .369) 

(.276, .322) 


, radians 

second 

(3.37, 11. «3) 


(5.27, 7.*1) 




deRrees 


(15.2, ftl.«i) 

(«I2.«I, 5*. 8) 


(36.7, 41.5) 

Parameters 
for Cro3!s- 
over Model, 

radians 
“*c* seconds 

3.2 

2.0 

1.9 

1.3 

1.5 

Based on 

Tg, seconds 

0.370 

0.533 

0.393 

0.572 

0.533 

Average 








1 

1 

i 
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where is the crossover frequency. Figure S shows a Bode diagram of 

jy)Y^; (jy) , where YpCjw) is the average of the estimated rider 
’ansfer functions from the second test day (from Table 1), and a 

iiode diagram of Yp( ju)Y^(j<a) for the 15 mph tests. In the 15 mph 
*ituation» it was*^necessary to use the average of Yc(j**>)» esti- 
mated from Equation (l)g because there was noticeably less agree- 
ment between the theoretical and the experimental as compared 
with the 30 roph tests. It can be seen that» with appropriate values 
of uiQ and the crossover model closely fits the experimental 
results near the crossover frequency. The values of 
were selected for each group of data are included in Table 1. 

As shown in Figure 3, |Yr(j«)| at 30 mph has a br i region 
of a *20 db/decade slope, with a steeper drop-off near 0 radians/ 
second. As speed decreases, this drop-off meves to lowjr fre- 
quencies. Thus, to fit the crossover model, increasing amounts of 
lead equalization are required as speed decreases. At 30 mph, 
however, lead is optional. 

Estimates of the remnant power spectra were made from the 
estimated rider impulse response functions. These calculated 
spectra, w^.ch, precisely speaking, were based on the estimates of 
Y (jw) befc.e they were fitted by Equations (2-4), were generally 
ot the form 




Constant 

25 ♦ (#2 


the ^orm that would be obtained if "white” noise were passed 
through a first order filte.’ having a break frequency of 5 radians/ 
second, which, in fact, was the filter applied to che data before 
digitizing. Thus, to at least about 10 radians/ second, the actual 
unfiltered remnant is apparently "white", as was assumed earlier* 

The estimated remnant consists of all the steering torque 
output of the rider tnat not related to the roll angle by the 
linear transfer function -Yp(ju). The remnant contains (1) rider 
output which is nonlinearly relater to tne roll angle, (2) path 
correction steering torques not linearly related to the roll angle 
(such as corrections needrd due to the real road not being per- 
fectly straight, obrtacle or uomp ^ *oidarce, etc.) , (3) miscellaneous 
steering torque#, voluntar/ or involuntary (including direct 
mechanical transmissioi of ro id shocks throi*gh the rider’s arms 
to the handlebars, although road shocks are high frequency 
phenomena and would not oe included in the estimated remnant), (4) any 
time variation in Yp(jw), and (5) errors in identifying Yp(jai). Since 
the rider is using torque control fox purposes other than roll 
stabilization, it was found thar the mean square remnant was a large 
proportion of the.anan square steering torque: from 30-95%. It is 



Figure 5. _ • 

Cx’ossover model approximations to Yp(Jto)Yg(j«) 

(rider A, 30 mph, second day average, Yp(3«) 

4 ?)), and (rider A, 

15 mph, first day average, Y^(Jm) • •75*6j«e ^ (^^ + D) 




not known how mucu oS this is dur to t-ror« in estimating Yp(j »). 
Probably the iuentif ication errors increased with the remnant, 
sin^e a remnant ended to .lake accurate identification of 

n*ore dif^icul.. 


4, CONCLUSIONS 

1. Uuiii)'’ normal op'^:ation o.i a paved road in good condition, 
the primary excitation to the i.an-m^ torcycte system is the rider's 
’’remnant**. This fact allows accurate identifier t ’on of the 
controlled elemer.^ CVc(j«j)) h; open-loop cross-spectral analysis, 

at least for the frequency range .n ^hich there is substantial 
amplitude of the records 1 eing ^^nalyzed. Agreement between the 
theoretical mu xperiir^ntal i better at 30 mph than at 

15 mph. 

2. The Wir^rove-E' ards metnod [5] can be used to identify 
Vp(ju>), the transfer functirn of "he ri^er for the task studied. 

^he experimental steering torqye power spectra were helpful in 
inte preting Bode diagrams of 

3. At 30 mph, :he . ide transfer function YpCjw) was found 

to be a consta.it ^nd t^me uelay, with the optional inclusion 

of lead equalization having a break frequency in the neighborhood 
of 5-10 "*adians/scc . ’*'ime 'delays i ere about 0.3 second for all 
the test., and gains were aoout iSC-350 Ib-in/radian. 

4. Only three experiments were performed at 15 mph, but the 
re^iults of these indi^cteu th«t a combination of rate ccutrol and 
lead equalization (break 'r'^quency, 5-10 radians /second) was 
required on the part of tne ridec. here the time delays were 

about 0.25-0.35 seco .d, and the gains were about 60-90 Ib-in/radians . 

5. The experimental results agree with the crossover model , 
v.hich has been used to thtoreticallv study the man-motorcycle 
svbtem [4]. 

6. Vest of the rider's steering torque output was remnant. 

The power spectrum cf n(t) was found to have about the same shape 
as that of *'white” noise passe i through the same filter that was 
appiied to the original analog t jcords , 
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APPENDIX A 


APPEN^^X B 


MEASUREMENT OF MOTORCVCLE PARAMETERS 


T'o calculate the theoretical Yc(jw) from the wotorcvcle 
equations of motion [1], it was necessary to ra<»ice a number o. 
neasurements of the test vehicle. 

In addition to geometric measurements, mass and mass distri- 
bution estimates were made for the front wheel, fork and handlebar 
assembly, and the total motorcycle with Rider A in place on the 
seat. fThe brace in Figure 2 was used.) The following quantities 
Kerr mea*;ured: 

1. Masses 

2. Center of gravity locations 

5. Roll moments of inertia (x-axi^J 

4. Yaw moments of inertia (t-axis) 

5. X2 products of inertia 

It was necessary to attach the front frame assembly and the entire 
vehicle to fixtures to measure nass distribution. Moments of inertia 
were measured with cab’ c* suspended torsional pendula, except for 
"he roll monent of inertia of the man-motorcycle combination, for which 
ineas jrements were made by swinging the assembly on knife edges. 

Addi t ioi ally , wheel polar moments of inertia were measured 
with a torsional ^/endulum arrangement. An effective polar moment 
of iner’^ia for tht engine anJ transmission was calculated from 
data supplied by the ma^.u acturer. 

Tire properties were obtained from the HSRI flat bed tire 
tester. Measurements included the dependence of tire lateral force 
and aligning r’^ments on slip and inclination aniiles, *^nd the tire 
relaxation length. 

Friction in the steeling head was estimated and found to be 
negligible, on the basis of analog computer studies not discussed 
here . 


IDENTIFICATION OF KNOWN THANSFER FUNCTIONS 

Artificial data were prepared by digital computer as follows: 
a subroutine generated a series of random numbers z(kh), k«l , 2, . . . 500 , 
h«0.1 second, where z(kh) was normally distributed with mean zero. 

The artificial remnant was calculated from 

n(kh^ = Xjnt(k-l)h] + 2(kh), k*l,2,...S00 

The theoretical autocorrelation function of n(kh) is [10] 

r„„(k'h) - aj k'-0,il,t2 

The ’’rider** and ’motorcycle*' transfer functions were taken to be 

Y_ (ju.) = 

Pi 

Yp^CjuO = 

and 

= h • 

The degree of bias in estimating Yp(j<«)) by the Wingrove-Edwards 
methou is increased by increasing rnn(Tp) =* rnn(0*^) [5], or, in 
thi^ case, by increasing Zero bias results if rnn(Tj®0 for 

T > 0.4. 

Figure 6 «hows the identification of Vp^ when * 0-3 
(^nn(0*^) Z ^P2 ^'l “ 

general, |Vp 2 l, which is mere biased, tends more toward 1/VcO*-) at 
the higher and lower fr^^quencies than |Ypj|. Further, 5(Yp2) is 

underestimated for u) ^ 2 radians/secord. If the bias is increased 
further, it has been found that both |tp| and HYp) tend to be 
underestimated. On the other hand, it is seen that |Yp^| and 

KYp ), which should not be strongly biased, tend to he slightly 
overestimated. 
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Figure 6, Identification of known "rider" transfer functions 
from artificial data by the time shifting method. 





